AD-A0A4  717 


UNCLASSIFIED 


NATIONAL  BUREAU  OF  STANDARDS  WASHINGTON  DC  QUANTUM  CH--ETC  F/6  20/s 
APPLICATION  OF  QUANTUM  CHEMISTRY  TO  ATMOSPHERIC  CHEMISTRY. <U> 

JUN  77  M KRAUSS.  F H MIES.  D B NEUMANN  AFOSR-ISSA-77-OOll 


»cn<;o-To-77-i  91  41 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  * STANDARDS -1963-r 


r— 1 

© 

< 

o 

<c 


/ 


U 


i 


AFdSRtTR-  7 7 - 12  14 

— - -J 

Semi-Annual  Repwrt.  FY1977,  | 


Application  of  Quantum  Chemistry 


to  Atmospheric  Chemistry  . 


Contract  No.  AF0SR-1SSA-77-OO11  j 

/ - 


Submitted  to: 

Air  Force  Office  of  Scientific  Research 
Bolling  Air  Force  Base 
Washington,  D.C.  20332 


n 3 


a_ 

O 

c_? 


Z3 


J M.  Krauss 
F.  H.  Mies 
D.  B.  Neumann 
P.  S.  Julienne 

Quantum  Chemistry  Group 
Physical  Chemistry  Division 
National  Bureau  of  Standards 
Washington,  D.C.  2023*4 


Approv  I 

distribution  u 


--■> : 


L;iii:ni  teu  > 


AIR  FORCE  OFFICE  OF  S0Ir-"'~' 
NOTICE  CF  T '■ 


■'  CH  (AFSC) 


l .1 


“it  on  Oi'iicor 


\ 


... 


— — 


REPORT  DOCUMENTATION  PAGE 


t.  REPOR"  NUMBER 


/ \2.  GOVT  ACCESSION  NO. 


AFOSH-TR-  7 7-12 


4 TITLE  (and  Subtitle) 

APPLICATION  OF  QUANTUM  CHEMISTRY  TO  ATMOSPHERIC 
CHEMISTRY 


KIM)  INSTKUCTIONS 
BEKOKE  COMPLETING  FORM 


3,  RECIPIENT'S  CATALOG  NUMBER 


5 TYPE  OF  REPORT  & PERIOD  COVERED 


Interim 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORfs; 

M Krauss  P S Julienne 

F H Mies 
D B Neumann 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Quantum  Chemistry  Group,  Physical  Chemistry  Div. 
National  Bureau  of  Standards 
Wash  DC  20234 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

AFOSR/NP 

Bolling  AFB,  Bldg. 410 
Wash  DC  20332 


14.  MONITORING  AGENCY  NAME  fit  ADDRESS^//  different  from  Controlling  Office) 


8.  CONTRACT  OR  GRANT  NUMBERfsj 


AFOSR-ISSA-77-OOl 1 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  4 WORK  UNIT  NUMBERS 


9767-01 

61102F 


12.  REPORT  DATE 

June  1977  / 


13.  NUMBER  OF  PAGES 

24 


15.  SECURITY  CLASS,  fof  thlt  roporf) 


16.  DISTRIBUTION  STATEMENT  (ol  this  Report) 


Unclassified 


1 5a.  DEC  L ASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 


Approved  for  public  release! 
distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  it  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  reverse  side  If  necessary  and  identify  by  block  number) 


20  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

The  coupled  Hartree-Fock(CHF) , method  was  applied  to  neutral  atoms  of  the  first 
row.  The  method  was  applied  with  and  without  the  inclusion  of  electron 
correlation.  When  electron  correlation  was  included  vis  the  coupled 
multiconf igurational  self-consistent  field,  (CMCSCF),  method,  which  is  a 
straightforward  extension  of  the  CHF  method  it  was  found  that  correlation 
effects  accound  for  less  than  10  percent  of  the  polarizability.  Further 
investigations  regarding  the  basis  set  via  including  mere  or  less  flexibility 
showed  that  extensions  of  basis  set  size  affected  the  polarizabilities  by  less 


DD  1 JAN  73  1473  EDITION  OF  1 NOV  65  IS  OBSOLETE 

SECURITY  CLASSIFICATION  OF  This  PAGE  rl*7,en  Pot*  Fn trr- 


9 U L'v 


m*  W*-\ 


X 


/ 


1%  Introduction 

X 

The  polarizabilities  of  ionic  states  of  atoms  are  quantities 
of  importance  in  the  determination  of  long  range  forces  between 

atoms  where  either  excited  states  are  present  which  are  ionic  in 

-j' 

character  or  which  exhibit  a large  degree  of  charge  transfer. 

The  electric  dipole  polarizability  represents  the  major  interaction 
of  a static  electric  field  with  an  atom  or  ion.  A considerable 
number  of  methods  have  been  proposed  and  applied  to  the  calculation 
of  atomic  polarizabilities.  These  range  from  the  solution  of  the 
perturbation  equations  describing  the  influence  of  the  field  on 
the  atoms  to  very  complete  variational  methods  with  large  numbers 
of  configurations. 

The  perturbation  equations  arising  when  atomic  exchange 

effects  are  neglected  have  been  solved  by  variational  techniques 

2 

for  small  atoms  by  Kirkwood,  Buckingham  and  others.  Dalgarno 
3 

and  Parkinson  achieved  moderate  accuracy  for  first  row  atoms  by 

solving  for  the  first  order  perturbed  orbitals  in  the  presence 

of  exchange.  More  recently  calculations  of  atomic  polarizabilities 

4-7 

have  been  performed  using  the  coupled  Hartree-Fock,  (CHF),  method, 

g_10 

the  coupled  Hartree-Fock  perturbation,  (CHFP),  method,  the 

11-13 

double  perturbation  theory,  and  the  many-body  perturbation 

theory. ^ ^ In  these  techniques  the  Hartree  or  Hartree-Fock  wave 
function  is  used  as  a zeroth-order  approximation  from  which  the 
first-order  perturbed  function  for  the  atom  in  the  field  is 

3 

obtained.  In  the  CHF  method  the  HF  equations  for  the  atom  are 
variationally  solved  in  the  presence  of  a finite  field.  In  the 


- 2 - 

CHFP  method  the  polarizability  is  expressed  as  a function  of  the 

3 

change  in  the  second  order  energy  of  the  atom,  while  in  the  CHF 
method  the  variationally  determined  wave  function  for  the  system 
in  the  presence  of  the  field  is  used  to  determine  the  induced 
dipole  moment  as  a function  of  the  field.  The  polarizability  may 
then  be  expressed  as  the  limit  of  the  ratio  of  the  induced  dipole 
moment  to  the  field  as  the  field  goes  to  zero.  If  sufficiently 
accurate  wave  functions  are  used  the  CHF  and  CHFP  methods  are 
equivalent . 

17  18 

In  previous  work  by  this  group  ’ the  CHF  method  was 
applied  to  neutral  atoms  of  the  first  row.  The  method  was 
applied  with  and  without  the  inclusion  of  electron  correlation. 

When  electron  correlation  was  included  via  the  coupled  multi- 
configurational  self-consistent  field,  CCMCSCF),  method,  which  is 
a straightforward  extension  of  the  CHF  method  it  was  found  that 

17  19 

correlation  effects  account  for  less  than  10%  of  the  polarizability.  * 

Further  investigations  regarding  the  basis  set  via  including  more 

or  less  flexibility  showed  that  extensions  of  basis  set  size 

18 

affected  the  polarizabilities  by  less  than  5%. 

In  the  present  work  a simple  ad  hoc  scheme  for  choosing 

basis  sets  of  ions  and  atoms  was  sought  which  when  used  with  a 

CHF  calculation,  without  the  inclusion  of  correlation  effects, 

would  give  polarizabilities  for  ions  to  10-20%  accuracy.  The 

accuracy  will  be  poorest  for  particular  atomic  states  where 

20 

correlation  effects  are  substantial. 

In  addition  to  the  work  on  the  polarizabilities,  we  have 

begun  to  prepare  the  recent  work  on  for  publication.  The 

abstract  for  the  paper  "The  + Transition  in  N " is  included 

u g 2 


II.  Review  of  the  Method 


In  the  previous  work  of  our  group  by  Billingsley,  Krauss  and 
17  18 

Stevens  * a CHF  method  was  used  in  which  the  atom  was  subjected 

to  static  field  due  to  a unit  charge  at  a relatively  large  distance 

10-20  a.u.  from  the  atom.  The  wave  function  in  the  presence  of 

the  charge  was  computed  variationally . From  the  dipole  moment 

induced  by  the  charge  the  polarizability  was  determined.  As  the 

field  is  not  uniform  the  dipole  is  induced  not  by  the  field  alone 

21 

but  also  by  gradients  and  higher  derivatives  of  the  field.  In 

the  present  work,  in  order  to  avoid  any  ambiguities  introduced  by 

the  nonuniformity  of  the  field  the  technique  of  Cohen  and  Roothaan 
4 5 

was  applied.  ’ In  their  application  of  the  CHF  technique  the 
atom  is  subjected  to  a uniform  field,  F.  The  total  Hamiltonian 
for  the  system  is: 

H = H + F h (1) 


N i 

where  h = -E  — = — and  H Is  the  Hamiltonian  for  the  unperturbed 

i«l  F 

ion.  The  field  used  is  directed  along  the  z-axis  of  the  system. 
This  destroys  the  spherical  symmetry  bringing  the  atom  to  C 
symmetry.  For  the  state  of  F this  produces  a ^E+  state.  It 
is  for  this  state  that  the  polarizability  is  determined.  The 
wave  function  can  be  written: 

ip  *=  ifi  + Fip 1 + F^ip  + • • • 
o J.  z 

where  F = |f|,  is  the  unperturbed  HF  wave  function  for  the 
N-electron  ion.  If  \p  were  the  exact  wave  function  of  the 


(2) 


unperturbed  system  then 


- A - 


(H  -E  )i|>  =0  (3) 

o o o 

and  ^ is  the  solution  of 

(H  -E  H.  + =0  (A) 

0 0 1 o 

In  such  a perturbation  approach  the  static  electric  dipole 
polarizability  becomes: 

a = -2<i|»  |h|i|^>  (5) 

In  the  present  CHF  method  the  wave  function  for  the  system  is 

determined  for  particular  values  of  the  field  along  the  z-axis 

and  one  can  calculate  the  dipole  moment  induced  by  the  field: 

V(F)  = <iji|h|i|i>  (6) 

The  polarizability  is  then  found  by  the  classical  definition: 

a « limly  (F)/FJ  (7) 

F-o 

when  higher  terms  are  included  then  the  induced  moment  may  be 
written  for  the  atom  or  ion  as: 

p(F)  = aF  + y F3  + * * * (8) 

where  y is  known  as  the  static  electric  dipole  hyperpolarizability. 

A second  method  of  determining  the  polarizability  is  from  the  change 
in  the  energy  cf  the  system  due  to  the  field.  The  energy  may  be 
written: 

E(F)  »=  «HHq  + Fh|^>  (9) 

- <*J>  + fr2{2<H-0|hU1>+<^1|Ho|*1>2<*o|HoU2>} 

2F4{  <ipQ  | h 1 4/3  >+<4^  | h | ^2>+<^i  I H0 1 ^3> 

+ ^2iH0iv+<^iHoiv} 

Terms  proportional  to  odd  powers  of  F are  absent  due  to  parity 
considerations.  Thus  by  using  the  relationship  for  the  exact 
first-order  wave  function  (A)  in  the  second  term  of  (9),  and 


J 


i 


and  omitting  the  probably  small  term  <1/;^  j Hq  | <^2> » one  °^tains: 
AE(F)  = E(F)-E(F=0) 

- F2<^o|h|^>+0(F4) 
which  when  coupled  with  (5)  gives 
ct(F)  = -2AE(F)/F2 

thus: 

a = lim{-2AE(F)/F^} 

F-*D 


(10) 

(ID 

(12) 


In  the  present  work  both  eqs.  (12)  and  (7)  were  applied.  It  was 

found  that  these  two  different  definitions  of  a do  not  give 

precisely  the  same  result.  This  is  basically  due  to  the  fact  that 

the  wave  functions  \li  , used  in  the  calculation  are  not  exact  as 

is  required  for  the  two  definitions  to  be  equivalent.  They  have 

5 22 

been  shown  also  to  be  equivalent  for  HF  wave  functions  by  Cohen  ’ 

with  the  further  proviso  that  Brillouins'  theorem  holds  for  ^ . 

Insofar  as  the  wave  functions,  <Jj  , used  in  this  work  are  not  exact 

o 

solutions  to  the  Hartree-Fock  equations  nor  does  Brillouins'  theorem 
hold  once  spherical  symmetry  is  broken,  it  is  not  surprising  that 
the  equivalence  between  the  two  definitions  is  found  not  to  hold 
exactly. 

In  addition  to  the  polarizability,  the  hyperpolarizability 
may  be  determined  by  fitting  the  induced  moment  (8)  or  from  the 
Initial  slope  of  y(F)/F.  In  a similar  manner  to  (5)  the  hyper- 
polarizability can  be  written  by  omitting  probably  small  terms: 

y m 12{<1^1 1 h 1 4»2>  + <iJ>o|h|tJ>3>}  (13) 


— » — MfcliaiHi  1 1 


I 
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III.  Computational  Methods 


The  calculations  reported  herein  were  performed  with  the 

23 


BISON-MC  system  of  programs  of  Wahl,  Das,  et  al.  The  BISON-MC 
programs  determine  variationally  the  wave  function  which  is 
expressed  as  a linear  combination  of  Slater  determinants.  In  the 
present  case  a single  determinant  was  used  as  correlation  effects 
were  ignored.  The  atomic  orbitals  of  the  Slater  determinant  are 
in  turn  expressed  as  linear  combinations  of  atomic  centered  Slater 
type  orbitals,  STO's.  The  wave  functions  were  generated  using 
the  Hamiltonian  Cl)  with  fields  varying  from  0.0020  to  0.0060  a.u. 


The  calculations  were  performed  in  symmetry.  Calculations 


with  the  field  set  to  zero  were  done.  These  gave  slightly  lower 


total  energies  than  the  calculations  done  in  spherical  symmetry. 


In  the  present  case  of  a S atom  this  lowering  was  quite  small. 


This  is  in  contrast  to  atomic  calculations  done  in  C^  symmetry 


17 


where  the  atom  was  in  a degenerate  state  where  larger  energy 
lowerings  occur  when  the  spherical  symmetry  is  broken. 


In  the  present  CHF  method  the  polarizability  is  calculated 
using  a basis  set  for  the  atom  consisting  of  two  parts.  The 
first  part  is  the  Hartree-Fock  basis  for  the  atom.  The  Hartree- 
Fock  basis  set  of  Clementi  for  F (^S)  was  used.  The  second  part 
of  the  basis  set  consists  of  the  polarization  functions.  Sitter 
and  Hurst  have  reported  rules  for  choosing  the  principal  quantum 
number  and  spherical  harmonic  portions  of  the  polarization  orbitals. 
There  remains  the  choice  of  the  exponents  of  these  polarization 


functions.  They  may  be  chosen  by  a method  described  by  Billingsley 
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and  Krauss  in  which  the  functional: 


*2  ■ <u:[1)ia°-s°l»'1,>  + 2<«ia,|h1|u«0,» 

is  minimized.  Here  is  the  Hartree-Fock  orbital  and  u^^  is 

the  polarizational  orbital. 

K - - i ”1 + W 


h = h h.  = F*r./F 

1=1  i’  l l' 


(14) 


(15) 


and  is  the  orbital  energy  of  the  i-th  orbital.  Though  this 
method  has  a firm  theoretical  base  it  is  difficult  to  apply  as 
the  potential  V is  a two  electron  potential;  hence  Billingsley 
and  Krauss  simplified  by  using  an  effective  potential: 

(o) 


W ■ 


* 4 * 

n (n  ^1) 


?i 


(16) 


r . 

l 
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due  to  Slater  and  Zener.  Still  in  order  to  apply  their  method  of 
choosing  exponents  for  the  polarization  functions  a further  simpli- 
fication was  made.  They  chose  to  represent  u^°^  and  u^^  by  single 
Slater  basis  functions. 

In  the  present  work  a similar  although  completely  ad  hoc 

method  was  used  to  choose  polarization  functions.  The  integral 

<u,(1)|h/|u<°>> 
i 1 i 1 x 

was  maximized,  and  u^°^  and  u(l)  were  represented  by  single  Slater 
orbitals.  This  gave  rise  to  a set  of  rules  to  be  used  in  combination 
with  those  of  Sitter  and  Hurst. ^ Thus  when  basis  functions  are  to 
be  chosen  to  allow  for  polarization  of  a Hartree-Fock  orbital 
whose  most  significant  function  is  Xn^°^  (C^°\  r)  choose  polarization 


JL  vlLi' 


functions  x^t’?»(C^>r)  by  the  rules: 
n X* 


(o). 


2s  (s 


(o). 


2p(C 


(o). 


with 

2Pa(1) 

« 

S(0)) 

and 

3Pa(1) 

= 

7s(o)/5) 

with 

2p(C(1) 

= 

5?  (o)/7) 

and 

3P(C(1) 

= 

.(0)) 

and 

4P«(1) 

= 

9C(o)/7) 

with 

isa(1) 

= 

3£(o)/7) 

and 

2s(^ 

= 

5S(0)/7) 

and 

3sU(1) 

= 

S(0)) 

and 

4s(s(1) 

= 

9^(o)/7) 

and 

3d(e(1) 

= 

;(0 

and 

4d(s(1) 

= 

9?(o)/7) 
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Should  ui(o)  have  more  than  one  significant  basis  function  then 
choose  polarization  functions  for  each  of  them.  Omit  any 
polarization  functions  which  strongly  overlap  the  space  of 
functions  used  to  describe  ip  to  prevent  over  completeness  of 
the  basis  set. 

In  using  the  Clementi  Hartree-Fock  basis  sets  for  an  atom*' 

it  frequently  occurs  that  there  are  three  or  more  functions  which 

have  large  atomic  orbital  coefficients.  In  such  a case  it  is 

necessary  to  pick  the  polarization  functions  using  the  above  rules 

but  based  on  the  exponents  of  a smaller  optimized  atomic  basis  set 

This  is  easily  done  (even  when  the  exponents  of  the  smaller  basis 

set  are  not  known  but  the  exponents  of  a very  small  basis  set  are 

25 

known)  by  an  extrapolation  technique  used  originally  by  Homback 


with  gaussians.  The  logarithm  of  the  exponents  is  plotted  versus 
the  size  of  the  basis  set  and  extrapolation  to  smaller  basis  sets 


is  performed.  In  figure  1 the  exponents  of  the  basis  functions 
for  the  2p  orbital  of  F are  plotted  for  several  basis  set  sizes. 

It  is  easily  seen  that  interpolation  or  extrapolation  should  be 
quite  accurate. 

Using  the  above  ad  hoc  rules  the  polarization  functions  for 
F were  chosen.  The  full  basis  set  is  given  in  Table  I as  Set  I. 
These  ad  hoc  rules  should  be  sufficient  to  compute  polarizabilities 
to  within  10  to  20%  accuracy  for  many  ions.  To  further  test  the 
method  of  choosing  the  basis  set  a second  larger  basis  set  was 
chosen  in  a manner  similar  to  that  of  Billingsley  ard  Krauss 
for  Carbon.^  This  basis  is  given  in  Table  I as  Set  II. 


m i maOmmoMk 
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IV • Results  for  F ( * S ) 

In  Table  II  the  change  in  the  energy,  AE(F),  relative  to 

E(F=0)  is  given  as  a function  of  the  field  for  basis  set  I. 

The  energy  for  the  spherically  symmetric  F ("S)  system  is  given 
26 

by  Clementi  as  -99.459363.  Using  basis  set  I in  the  absence 

of  a field  we  obtain  a lowering  of  the  energy  by  only  5.4822x10  ^ a.u. 

This  small  change  is  due  to  the  inclusion  of  the  2s  polarization 

function.  In  Table  II  the  induced  dipole  moment  as  a function  of 

the  field  is  also  presented.  Due  to  the  use  of  a single  component 

of  the  field  F rather  than  three  equal  field  components  F = F 
z x y 

= F a slight  quadrupole  moment  is  induced  into  the  system.  This 
z 

1 2 2 

moment  - | I (3z^-r^) | i|>>  is  given  in  Table  II  as  a function  of 

i 

the  field  as  is  also  the  relative  value,  AX,  of  the  expectation 
N 2 2 

value  of  Z (xrty\).  The  polarizability  as  a function  of  the  field 
is  given  by  both  the  expressions  (7)  and  (11).  These  values  are 
plotted  in  Fig.  2.  From  Table  II  it  can  also  be  seen  that  the 
small  quadrupole  moment  induced  in  the  ion  is  directly  proportional 
to  F with  a slight  higher  order  effect  present.  Thus  from  Table  II 
it  can  be  seen  that  both  AX  and  the  magnitude  of  the  quadrupole  are 
increasing  with  the  field,  but  the  overall  changes  are  very  small. 

From  the  plot  of  the  polarizability  as  a function  of  the  field  in 
Fig.  2,  it  can  be  seen  that  basis  set  I gives  a static  dipole  polarizability 
of  9.398  a.u.  using  the  induced  dipole  moment  expression  and  a slightly 
larger  value  9.402  a.u.  using  the  energy  expression. 


v»*>  fa*  < 
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In  Table  III  the  results  are  given  for  the  calculation  using 
basis  set  II.  The  polarizability  via  expressions  (7)  and  (11) 
is  plotted  for  basis  set  II  in  Fig.  2.  The  dipole  moment 
expression,  (7),  gives  a polarizability  of  8.804  a.u.  while 
the  energy  expression,  (11),  gives  a polarizability  slightly 
larger  8.810  a.u. 

The  plots  of  polarizability  versus  the  square  of  the  field, 

as  shown  in  Fig.  2,  show  that  the  limiting  value  of  a is  only 

slightly  different  when  expressions  (7)  or  (11)  are  used.  Second, 

the  plot  shows  especially  for  basis  set  II  a certain  amount  of 

scatter  when  the  induced  dipole  moment  is  used  via  expression  (7) 

to  determine  the  polarizability.  There  seems  to  be  considerably 

less  scatter  when  the  second  order  energy  is  used  via  expression 

(11)  for  the  polarizability.  The  explanation  for  this  is  that  the 

self-consistent-field  algorithm  used  to  determine  the  wave  function 

gives  a wave  function  which  is  self-consistent  or  precise  to  an 

extent  e,  wThile  the  total  energy  determined  in  this  method  has  an 

accuracy  which  goes  as  the  square  of  the  accuracy  of  the  vectors 
2 

i.e.  as  e . The  dipole  moment  accuracy  follows  that  of  the  vectors. 
These  statements  can  be  shown  to  follow  rigorously  as  a consequence 
of  the  variational  method  used  and  as  a consequence  of  Brillouins* 
theorem.  At  a latter  stage  in  the  present  investigations  we  intend 
to  examine  this  problem  of  the  accuracy  of  the  method  and  the 
applicability  of  Brillouins' theorem,  especially  for  open  shell 
atomic  systems. 
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These  results  are  to  be  compared  with  the  results  of  other 

calculations  for  the  polarizability  of  F (^S)  as  given  in  Table  IV. 

The  comparison  shows  that  the  polarizabilities  calculated  by  the 

straightforward  CHF  procedure  via  a relatively  ad  hoc  scheme  of 

choosing  the  basis  set  agrees  fairly  well  with  the  results  of  other 

methods  at  the  same  level  of  accuracy. 

This  procedure  is  now  being  applied  to  positive  and  negative 

ion  states  with  s or  p orbital  open  shells  as  well  as  closed 

shell  systems.  At  the  present  time  there  is  very  little  work  on 
27 

open  shell  systems  of  ions. 
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V.  Abstracts  of  Publications 

The  5I1  -^Z+  Transition  in  N„,  submitted  to  Mol.  Phys. 
u g 2’ 

Accurate  electronic  energy  curves  and  wave  functions  of  1 n 
and  1^Z+  states  of  N„  have  been  calculated  using  the  multi-configuration 

g 2 

self-consistent-field  (MC-SCF)  method.  The  state  was  predicted 

by  Mulliken,  who  suggested  it  was  important  in  afterglow  emissions. 

The  calculated  R and  D of  the  ^IT  state  is  in  reasonable  accord 
e e u 

with  estimates  of  Carroll  and  Mulliken  based  on  predissociative 
3 

behavior  of  the  C IT  state. 


u 

The  vertical  transition  energy,  is  calculated  to  be 

about  1.8-1. 9 eV.  The  bound  and  resonance  vibrational  structure 

of  the  state  are  such  that  the  transition  would  be  observable 
g 

as  bound-like.  The  calculated  transition  probability  is  zero 
asymptotically  but  increases  rapidly  to  shorter  distances  as  the 
atoms  overlap.  At  the  vertical  transition  distance  the  transition 
probability  is  about  2*10^s  " or  a lifetime  of  5usec  for  the  ground 
vibrational  level. 


- 
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(a)  For  functions  with  quantum  number  i 

greater  than  zero  both  a and  it  components 
were  included. 


(b)  Hartree-Fock  atomic  basis  taken  from 
Ref.  26. 
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SCPT:  Self  consistent  perturbation  theory 
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Figure  1.  Plot  of  exponents,  5,  of  STO's  used  to  represent 
the  2p  orbital  of  F~(^S)  as  a function  of  basis 
set  size.  Basis  functions  from  Ref.  23. 

Figure  2.  Plot  of  o(F),  the  polarizability  as  a function  of 

the  square  of  the  applied  field.  Plot  A and  B apply 
to  basis  set  I,  using  expressions  (7)  and  (11), 
respectively.  P’ot  C and  D apply  to  basis  set  II, 
using  expressions  (7)  and  (11),  respectively.  Units 
are  atomic  units;  conversions  are  given  in  Table  II. 
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